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Abstract 

(N 

We present a novel approach for computing the surface roughness-limited thermal conductiv- 
ity of silicon nanowires with diameter D < 100 nm. A frequency-dependent phonon scattering 
rate is computed from perturbation theory and related to a description of the surface through the 
root-mean-square roughness height A and autocovariance length L. Using a full-phonon dispersion 

g 

i relation, we find a quadratic dependence of thermal conductivity on diameter and roughness as 
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(D/A) 2 . Computed results show excellent agreement with experimental data for a wide diame- 
i ter and temperature range (25-350 K), and successfully predict the extraordinarily low thermal 
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conductivity of 2 W m K at room temperature in rough-etched 50 nm silicon nanowires. 
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Silicon nanowires (NW) have drawn much attention for their potential applications in field 
effect transistors [I], interconnects [2], thermoelectrics |31 IH El El E], and heterostructures 
[8]. Given their high surface-to-volume aspect ratio, the most prominent size effect limiting 
transport originates from electron or phonon surface scattering. In particular, being able to 
independently control electrical and thermal conductivity of such nanostructures through 
geometry, strain, or doping is extremely appealing for novel applications such as thermo- 
electrics and energy transport. Earlier data shows that reducing the NW diameter below 
100 nm leads to a drastic reduction in their electronic and thermal conductivity [2J El U\- 
More puzzling are the recent experimental results of Hochbaum and Boukai |3J 0] which 
show that intentionally etched rough edges reduce the thermal conductivity of crystalline 
silicon NW by a factor of about 100, to nearly the value of amorphous silicon. Several efforts 
have been previously made toward an accurate understanding of phonon-surface scattering 
[SI El IB], however, no studies account for such experimental observations in very rough wires 
with diameter below 50 nm. Similarly, no model provides guidance on how the thermal 
conductivity of such NW scales with their surface roughness. 

In this Letter we introduce a comprehensive approach to phonon-surface scattering in thin 
NW based on a perturbative treatment of interface roughness. We derive a matrix element 
for phonon-surface scattering which is directly related to a parametric description of surface 
roughness. Based on this approach, the effects of temperature change and surface quality 
on the phonon-surface scattering rate are evaluated in silicon NW of diameters below 100 
nm. Resulting theoretical predictions of NW thermal conductivity shows excellent agreement 
with experimental values below 50 nm, where the effect of surface roughness is the strongest. 

When the characteristic dimensions of asperities at a rough surface come to the order of 
the phonon or electron wavelength (5-30 A), it is expected that the surface scattering rate 
will be altered to reflect the effect of the interface roughness. Such mechanisms are crucial 
to the understanding of electron transport in transistor inversion layers, where accurate 
models have been developed based on perturbation theory [10] ■ While such formalism 
exists in the case of bulk phonon transport [TT] . current models of phonon-surface scattering 
in NW are based on simplified assumptions, most of them using the probability of diffuse 
scattering as a fitting parameter [2j [6]. Yet, the latter probability can be directly related 
to physical properties of the interface, which, among other options, may be experimentally 
observed by means of Transmission Electron Microscopy (TEM). Besides, it seems relevant 



2 



that the effect of rough surfaces should be stronger in thin NW, and vary with the frequency 
of incident phonons. 

In a thin nanowire, variations of the confinement width perpendicular to the propagation 
direction influence phonon transport by perturbing the Hamiltonian of the system (Fig. 
[T] (a) inset). It is assumed that boundary scattering is mainly an elastic process and no 
phonons are emitted into the surrounding environment. This condition reflects the case 
where NW are wrapped in a medium of considerably different thermal conductivity, as it is 
for Si NW in Si0 2 or vacuum. We model phonon transport in such Si NW of diameter below 
115 nm, comparable to experimentally available data. We introduce a new type of phonon 
scattering originating from the roughness of the NW surface. In essence, this scattering 
mechanism accounts for the fact that phonons "see" a rough NW as a series of constrictions 
along their propagation direction. In order to accurately model this effect at the nanometer 
scale, perturbation theory is used to derive the transition probability per unit time of an 
incident phonon of momentum k and energy fuv to a new state of momentum k' and energy 
ho' due to the perturbed Hamiltonian H' [12]: 

Under a sufficiently long time t in comparison to the energy relaxation time, the time deriva- 
tive reduces to the Dirac delta function. The interface roughness is considered as a space 
varying dilation A(r) of the wire. This alters frequencies in a plane perpendicular to the 
propagation direction in such a way that u;(k) = u; (k)[l — 7 A] where 7 is a fitting constant 
determined from thermal expansion of the material, and cuo(k) is the phonon dispersion of 
the unperturbed Hamiltonian. Following the derivation of Klemens [11], the matrix element 
for a perturbation due to a space varying dilation is 

l(k|^|k')| 2 = ^' 2 ((n) + l)A(k - k') (2) 

where V Q \ is the volume of the device, and A(q) is the Fourier transform of the spatial 
perturbation, equal to A(q) = f A(r)e* q ' r dr. Additionally, the occupation number is given 
by the Bose-Einstein distribution and includes the temperature dependence of the scattering 
process (n) = ( e hw / kBT — l) . 

As shown by Goodnick et al. [13], the autocovariance function of Si surface roughness is 
"roughly" fit by a Gaussian function, which, by the Wiener-Khinchin theorem yields a power 
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spectrum of 

A(q) = 7rA 2 LV> 2L2 / 4 (3) 

where A is the root-mean-square (rms) value of the roughness fluctuations and L is the 
autocovariance length, which is related to the mean distance between roughness peaks at 
the Si-Si0 2 interface (see Fig. [T](a) inset). In practice, these process-dependent parameters 
are experimentally set by the quality of the surface. Due to the u' 2 term in equation [2j 
low frequency phonons see little contribution from the surface perturbation. On the other 
hand, the power spectrum of equation [3] favors scattering processes of the specular type. 
Hence, one can expect that there is a frequency range over which phonons experience a 
higher contribution from surface roughness scattering. The phonon scattering rate from a 
branch i to a branch j is given by 

r-j (E) = J P h3 (k,k')dk' (4) 

The volume integral over k'-space can be reduced to a surface integral as shown in [14J: 

T - m ^j^f «L S (5) 



where Ni(E) is the phonon density of states in the i t h branch, and E'(k') goes along the j t h 
branch. The total scattering rate Ti(E) starting in branch i is the sum over all branches j 
of the Tij(E). 

A Gilat-Raubenheimer (GR) scheme [15J is used to compute surface integrals, which 
constitutes the optimal trade-off between accuracy and computational efficiency. In order 
to carefully account for frequency dependence, a full phonon dispersion is used, which is 
obtained from an adiabatic bond charge model and tabulated for look-up [El H7]. The 
GR method is also similarly applied to compute the phonon density of states based on the 
dispersion relation mentioned above. With this respect, the GR scheme divides the first 
Brillouin zone in a lattice of 40x40x40 cubes, which achieves sufficient accuracy in the 
10-350 K range p2J. Due to Bose-Einstein statistics, the accuracy of the scheme decays at 
lower temperatures, where finer grids or analytical derivation are required. Modeled silicon 
NW have a square cross section with area equivalent to a circular cross section of diameter 
D } while the NW length is arbitrarily fixed to 2 /xm. Transitions among all acoustic and 
optical branches are considered. 
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The scattering rates are first computed for NW of equivalent diameter D = 115 nm at T = 
300 K. For a fixed correlation length L = 6 nm, phonon lifetime is calculated for increasing 
A in Fig. [ijb). Although no reliable data on the roughness of small wires is available yet, 
various theoretical and experimental studies have reported roughness rms ranging from 3 A 
to 5 nm in the case of extremely rough NW pH HE]. In this study we consider effective values 
for A in this range, and first use the 3 A value for "smooth" NW, while using an average of 
3 nm for "rough" NW. Under the conditions cited above, the average phonon lifetime due 
to surface roughness alone is approximately 15 ps, and decreases with higher roughness rms 
values as A~ 2 . As deduced from equation [3j a long correlation length L favors scattering 
processes close to the specular type. While in the strong roughness limit where A/L > 1 the 
Gaussian approximation may be put at fault, the effect of the L 2 term tends to average out 
the contribution of L. We noticed only little deviation of the predicted thermal conductivity 
in the strong roughness limit, and consistently used a value of L = 6 nm, which is estimated 
from the TEM images of Hochbaum et al. |3] and provides a best fit in our case. 

In our approach, an additional thermal variation of the phonon-surface scattering rate 
appears from Bose-Einstein statistics, where temperature delimits the occupation of each 
frequency range, thus retrieving the effect of occupation of lower energy branches at low 
temperature. Subsequently, it is possible to determine the thermal conductivity for NW of 
different cross sections [19]. The contribution to thermal conductivity of branch i is 

Ki {T) = 1 jENiE^vmrnE) dE (6) 

where t>j is the velocity of sound which is dependent on the direction of propagation, here 
assumed to be in the <001> direction. In order to reproduce the measured physical behav- 
ior of the NW in the 10-350 K temperature range, Umklapp, normal, impurity, boundary, 
and surface roughness scattering mechanisms have been considered in the derivation of the 
branch-specific scattering time, as summarized in Table [T| Umklapp scattering in transverse 
acoustic branches is efficiently described with the law derived in Ref. [20], which has shown 
good agreement with Si NW experiments. Additionally, in the temperature range consid- 
ered, normal scattering in longitudinal acoustic branches is accounted for according to the 
derivation of Holland [21] • A Cu 4 law is often used in the literature for impurity scattering, 
and we found that a constant C = 8xl0 -45 s/K 4 fits experimental data with more accuracy 
than the value of 1.05 xl0~ 44 s/K 4 resulting from the derivation of Ref [21 J . In the case of 
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NW etched from a 10 Q-cm wafer |3], the low doping concentration has almost no impact 
on the value of C. Besides, the boundary scattering rate depends on the sound velocity in 
a given branch, and is consequently frequency dependent. An analytical expression is used 
for the sound velocity in acoustic branches, which is derived from [22] for its good fit with 
bulk Si data in the <001> direction. 

Fig. [2] (a) and [2] (b) compare thermal conductivity computed from the model presented 
above vs. experimental data from ['5| and [3], for NW ranging from 115 nm to 22 nm diam- 
eter. As nothing guarantees that all NW have similar roughness parameters, we considered 
smooth NW with 1 A < A < 3 A, and rough NW with 3 nm < A < 3.25 nm. A good fit is 
found for smooth NW grown by vapor-liquid-solid (VLS) mechanism with diameter above 
37 nm. Similarly, the model reproduces the drastic decrease in thermal conductivity for 
rough electroless-etching (EE) NW presented in |3]. Higher discrepancy is found for the 22 
nm smooth VLS NW, for which the sensitivity to surface roughness is expected to be higher. 
Besides, the perturbative approach remains valid as long as perturbations remain small in 
comparison to the total phonon energy. For low temperature phonons and nanowires of di- 
ameter below 20 nm, explicit quantum treatment may be required. Since surface roughness 
scattering has little impact on low energy phonons, additional low temperature discrepancy 
is attributed to impurity and classical boundary scattering. 

It is important to point out that our approach based on perturbation theory introduces 
a dependence in (D/A) 2 of the NW thermal conductivity, in contrast to the typical linear 
scaling with D used in previous descriptions. In this scope, figures [3] (a) and [3] (b) show 
how the thermal conductivity is lowered by concurrent effects of small diameter and rough 
surfaces. In particular, there exists a critical diameter below which the roughness-limited 
thermal conductivity of the NW noticeably deviates from the classical linear approximation. 
The fact that this critical diameter increases with higher roughness rms A, in turn, supports 
the expectation that heat conduction at small NW scales is strongly limited by their surface 
roughness. This assumption is further justified in Fig. |4[ where it is observed that the total 
contribution of surface roughness to limiting the thermal conductivity is increased about 
7 times from the 115 nm to the 22 nm case. Finally, the temperature dependence that 
arises from the roughness scattering rate is also retrieved. Thus, at low temperatures, low 
frequency phonons see little effect from the perturbations at the silicon interface. 

In summary, we have demonstrated a remarkably strong effect of surface roughness on the 
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thermal conductivity of thin silicon nanowires. Based on a full phonon dispersion relation 
we introduced a novel frequency-dependent model of boundary scattering for phonons. The 
resulting simulated thermal conductivity of NW of diameter D < 115 nm shows excellent 
agreement with recent experimental work. In particular, at low NW diameters, we predict 
a strong deviation of the roughness-limited thermal conductivity from the linear diameter 
dependence (~ D) to a scaling as (D/A) 2 , where A is the rms surface roughness. The 
approach presented here can be generally extended to NW of various materials, sizes, and 
direction of thermal propagation using the same full-phonon dispersion relation with the 
inclusion of the surface-scattering rate in a Monte Carlo simulation. 

This work was supported in part by the DARPA through the IMPACT Center for Ad- 
vancement of MEMS/NEMS VLSI grant n. HR001 1-06-1-0046 (PM), the DOE Compu- 
tational Science Graduate Fellowship Program of the OSNNSA under contract DE-FG02- 
97ER25308 (ZA), and the Nanoelectronics Research Initiative (NRI) SWAN center (EP). 
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Mechanism 


Analytical Model 


Constants 


Phonon-Phonon: 






Longitudinal [21] 


B L T 3 u; 2 


5 L =2xlO" 24 s/K 3 


Transverse [2"U] 


B T Te-^/ T iu 2 


B r =l-73 xl0~ 19 s/K 






0i = 137.3 K 


Impurity 121 


Cuj 4 


C=8xl0" 45 s/K 4 


Boundary 


v iD^ 1 a/1 + AujCi/vi 


^ L =9.01xl0 3 m/s 


nana 




v T =5.23xl0 3 m/s 






c L =-2xl(T 7 m/s 2 






c T =-2.26xlO" 7 m/s 2 


Surf. Roughness 


cf. Eq. H 


A= 3 - 50 A 






L = 60 A 



TABLE I: Summary of the scattering processes considered in the derivation of the thermal conduc- 
tivity. 
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FIG. 1: a) Phonon density of states for silicon computed from a full phonon dispersion relation, b) 
Phonon relaxation time due to interface roughness for rms A = 1 nm, 2 nm, and 3 nm, T = 300 
K, and a cross section equivalent to a circle of 115 nm diameter. L is fixed to 6 nm. 
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FIG. 2: (a) Thermal conductivity of smooth VLS Si NW. Shaded areas are theoretical predictions 
with roughness r.m.s. A = 1 - 3 A, blue squares are taken from [5], (b) Thermal conductivity of 
rough EE Si NW (A = 3 - 3.25 nm). Squares are taken from [3J. Simulation and experimental data 
are compared at similar cross sections. L = 6 nm. 
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FIG. 3: a) Predicted effect of roughness rms on the thermal conductivity of 115 nm and 56 nm 
NW at T = 300 K, L = 6 nm. Simulation and experimental data are compared at similar cross 
sections, (b) Effect of NW equivalent circular diameter on thermal conductivity [L = 6 nm). 
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FIG. 4: Proportional contribution of the surface roughness scattering term to the thermal conduc- 
tivity of VLS NW (A = 3 A, L = 6 nm). The reference value K* is the thermal conductivity of a 
perfectly smooth wire (A = nm). 
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